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Abstract 

5D models with one 3D brane and one infinite extra dimension are studied. Matter is 
confined to the brane, gravity extends to the bulk. Models with positive and negative tension 
of the brane are studied. Cosmological solutions on the brane are obtained by solving the 
generalized Friedmann equation. As the input in cosmological solutions we use the present-time 
observational cosmological parameters. We find constraints on dimensionless combinations of 
scales of 5D models which follow from the requirement that with cosmological solutions we 
reproduce the data on production of 4 He in primordial nucleosynthesis. 

1 Introduction 

Models with extra dimensions naturally appear in microscopic theories unifying gravity with other 
interactions. Viewed as models of the real world they must be tested as cosmological models. In 
this work we consider some aspects of cosmological solutions in five-dimensional with one 3D-brane 
embedded in the bulk with one infinite extra dimension [1] . Matter is confined to the brane, and 5D 
gravity extends to the bulk. 

There are no strong constraints on scales of the models with extra dimensions. In the braneworld 
DGP models [2] cosmology was discussed in a number of papers (the review article [3], also [4, 5, 
6, 7, 8, 9]). In these studies the fundamental 5D scale varied in the range 10~ 12 GeV — lGeV. The 
lower bound is provided by the tests of the Newton law, the upper bound follows from cosmology. 
Because the action of the DGP model does not contain cosmological constants in the bulk and on the 
brane, the metrics obtained by solving the Einstein equations contain no warp factors. The metrics 
considered in these models were taken static and either asymptotically flat [6, 7, 8], or deSitter [9]. 

In this note we consider models which actions contain both bulk and brane cosmological constants. 
In the Gaussian normal frame metrics in the bulk are warped solutions of the Einstein equations, 
satisfying Israel junction conditions on the brane with matter [10, 11, 12, 13, 14]. The warp factor 
is taken in the range 10~ 12 GeV — 10 3 GeV. 

In models with extra dimension cosmological evolution on the brane is described by solutions of 
the non-standard Friedmann equation [11, 12, 13, 14, 6]. We do not make a fit of parameters of 
models using the observational cosmological data, but taking as the input the set of present-time 
observational cosmological parameters (Hubble parameter, deceleration parameter and fractions of 
cold matter and radiation in the total energy density) we look for constraints on parameters of the 
models which follow from the requirement that the models yield abundance of 4 ife produced in 
primordial nucleosynthesis within the experimental bounds. 
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We consider models with positive and negative tensions of the brane. The latter case is suggested 
by models with two branes, which allow for a possibility to reconsider the hierarchy problem. In 
these models the visible brane should have negative tension [15]. 

Calculating production of 4 He in the Big Bang nucleosynthesis (BBN) we find constraints on 
dimensionless combinations of scales of the models. In particular, we find a relation 

where \x and M are the scale of the warp factor of the metric and 5D fundamental scale, and r c M 3 
is the coupling at the 4D Einstein term of the action. 



2 One-brane model 

We study one-brane models embedded in 5D bulk with the action 

W^(S +a W^(§-*)-/^ l "" (i) 

where k 2 = 8n/M 3 , k\ = 8n/r c M 3 . L m is the Lagrangian of matter on the brane. M is the 
gravitational scale of the 5D gravity, parameter r c defines the strength of gravity in the 4D term. 

It is assumed that matter is confined to the spatially flat brane and gravity is 5-dimensional. To 
have homogeneous cosmology on the brane the energy-momentum tensor of matter on the brane is 
taken in a phenomenological form 

T£ = diag{-p,p,p,p}, (2) 

where pit) and p(t) are the sums of energy densities and pressures of cold matter and radiation. 
For the following it is convenient to introduce the normalized expressions for bulk cosmological 
constant, energy density, pressure and cosmological constant on the brane which all have the same 
dimensionality GeV 

I K 2 A K 2 p K 2 a n 2 p . 

The Einstein equations are solved in a class of metrics of the form 

dsl = -n 2 (y, t)dt 2 + a 2 (y, t)-f ij dx i dx j + dy 2 . (4) 

The brane is located at the fixed position y = 0. Using the freedom in parametrization of time, the 
component of the metric n(y,t) is normalized so that n(0,t) = 1. 

We study cosmologies on the brane by solving the non-standard Friedmann equation [11, 12, 14, 6] 
which follows from the system of the Einstein equations and junction (Israel) conditions on the brane 
(in the following we assume the symmetry y — > —y) 

a'(0,*) V- m , r c d 2 (0,t) 

n'(0,t) ^ r c ( d 2 (0,t) o(M) \ 



where pi, Pi i = r,m are densities and pressures of the cold matter and radiation. Equation for the 
scale factor a(0,t) can be obtained either in the form with the second-order derivatives of the scale 
factor [11], 

or, starting from the partially integrated system of the Einstein equations, in the form with the 
first-order derivatives [12, 13, 14, 6] 



d 2 (0,t) 



r c d*(0,t) \ 2 a 4 ((Mo) 



(7) 



a 2 (0,t) 

where the last term is interpreted as dark radiation [1, 16]. Introducing notations 

z = a(0,t )/a(0,t)-l, H(t)=a(0,t)/a(0,t), C{z + l) 4 = pp w (z), 
where to is the present time, Eq.(7) can be written as 

/ r - c H 2 \ 2 

(1 + ar c )H 2 = a 2 - p 2 + 2a(p m + p r ) + f p m + p r - J + pp w (8) 

Energy densities of cold matter and radiation are 

Pm = Pmo(l + zf, pr = PrOr(z)(l + Z) 4 ,. (9) 

where r(z) is a slow function of z which counts the number of relativistic degrees of freedom, r(0) = 1. 
If expressed as a function of temperature of the Universe, p r (T) = n 2 / 30 g*(T)T 4 [17]. This makes 
possible to connect r(z) and g*(T). 

2.1 Model with positive tension of the brane 

We define the period of "late cosmology" as a period in which in the Friedmann Eq.(8) the terms 
linear in matter (radiation) energy density are dominant 

ap(z) > p 2 (z), ap(z) > r c H 2 p(z), ap(z) > (r c H 2 ) 2 . (10) 

Friedmann Eq. (8) can be approximately written as 

(1 + r c a)H 2 ~a 2 -p 2 + 2a(p m (z) + p r {z)) + pp w (z). (11) 

Using (11) it can be shown that the second and the third conditions (10) follow from the first 
condition a > p(z). 

From Eqs. (6) and (7) taken at present time we express a 2 — p? and C . Introducing dimensionless 
combination 

P 2 = (12) 
and neglecting the terms of the second order in the present-time energy densities, we have 

a 2 -p 2 ^^[(l + r c a)(l-q )p 2 -l] : (13) 



and 



2p 2 



-3-A^ + (l + r c a)(l + q )p 2 



(14) 



where go = —d(to)/a(to)HQ is the present-time deceleration parameter. 

Substituting in (11) expressions (13) and (14), we obtain Friedmann equation in the period of 
late cosmology 



H, 



H 2 ^^>[l- qo + (l + qo )(z + iy 



(15) 



p 2 (l + ar c ) 
The rhs of (15) is positive if 



-1 + 4(s + If - ( 3 + 4^p(r(z) - 1) ) 1) 



p 2 (l + (7r c )(l + go) >3. 



We assume that parameter r c is constrained so that 

H 2 

Hr c -± < 1. 
/i 2 



(16) 



(17) 



For the smallest fi ~ 10 12 GeU we obtain the estimate /xr c < 10 60 . Below it will also be shown that 
p 2 {\ + or c )(l + go) — 3 <C 1, and hence 



p 2 ^ 



(l + <rr c )(l + g )' 
Using the constraints (16) and (17) in equation (13), we obtain that 



(18) 



a = ji + O ( jir c 



(19) 



With this accuracy we substitute /i for a. At large z the first condition of late cosmology (10) is 

/ \ 4 / \ 4 

// > p r0 r(z)2; = p m0 — r(z)z . 



Using (12) and (18) this condition can be written as 

z r(z) < 



O-iJ, 2 fi r # 2 (l +/xr c )' 
In the model with r c = 0, in the radiation-dominated period, Friedmann equation is 



2p 2 



n 



H 2 = ^[-3 + 4^(r(z) - 1) + (1 + q )p 2 z 4 + 



ij2 n \ 2 

— — r(z)z 
p 2 pil m 



where we have kept both the linear and quadratic terms. The linear term is dominant, if 

1 > ^la r r 2 (z)z\ 



(20) 



(21) 
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which is the same as condition of late cosmology (20) 

Because the present-time total energy density of the Universe is approximately equal to the critical 
energy density p c = 3ifoM^/87r, we can express p 2 as 

2 H 2 2M 3 

P — = 5 

pVt m p c tomflM^ ' 

where p c = K 2 p c /6. From the relation p 2 (pr c + 1) ~ 3/(1 + g ) we obtain 

/fjg _ 2(1 + g ) 
M 3 ~ 3^ m 

Here the combination 3Q m /2(l + g ) is of order unity 1 . 



(pr c + l). (23) 



3 Model with negative tension of the brane 

In the model with negative tension of the brane we set in (6) and (7) a — — \a\. The period of late 
cosmology is defined by conditions 

\a\p(z) > p 2 (z), \a\p(z) > r c H 2 P (z), \a\p(z) > [r c H 2 ) 2 

ensuring that in the Friedmann equation the terms linear in energy densities are dominant. Taking 
Eqs. (6) and (7) at present time and omitting the terms of the second order in energy densities, we 
have 

H 2 (l - r c \a\)(l - g ) ^ 2(a 2 - p 2 ) - \a\p m0 , (24) 
H 2 (l - r c \a\) - -V 2 + ° 2 - 2M(p m0 + Pro) + C. (25) 

From these relations we express a 2 — p 2 and C 

a 2 -p 2 ^^l(l + (l- g )(l - \a\r c )p 2 ) , (26) 



2 



2p 2 



2p 

3 + 4^ r - + (l-|a|r c )(l + g )p 2 

1 L rr) _ 



(27) 



In the model with r c = 0, in the radiation-dominated period, keeping both the terms linear and 
quadratic in the radiation energy density we obtain the Friedmann equation as 



tf- (3 + (1 + + , (28) 

where we neglected the term (r(z) — l)Vt r /Vt m as compared to 3 + (1 + qo)p 2 - 



1 With q ~ -1 + 3f2 m /2 [18, 20, 19], one obtains 3f2 m /2(l + q ) ~ 1. 



In the model with r c 7^ 0, neglecting quadratic terms and setting r(z) = 1, the Friedmann 
equation in the period of late cosmology is 



2 



1 - g ) + (1 + «b)(l + zf + p2(1 _ 1 Hre) (! " 4(« + I) 3 + 3(z + l) 4 ) 



(29) 



In the case 1 — |<r|r c > 0, for |g | < 1 the rhs of (29) is a positive increasing function for all z > 0, 
and there are no constraints on p 2 . 

In the case \cr\r c — 1 > the Friedmann equation in the period of late cosmology can be written 



as 



H, 



I! - t£ 

2 



1 



1 - 9o) + (1 + <Zo)(l + zf + -^th tt ("I + ±(z + I) 3 - 3(z + l) 4 ) 

p z (\a\r r . - 1) 



(30) 



Equation (30) has the same functional form as that in the model with positive tension. The rhs of 
(30) is positive for 

p 2 ( /U r c -l)(l + g )>3. (31) 
Constraining \ir c so that 1 3> {jir c )Hl / ji 2 , we can show that 



\a\ =ii + o(^rJ^ . 



4 Primordial nucleosynthesis 

4.1 Model with positive tension of the brane 

In this section, comparing predictions of BBN in the standard and non-standard cosmologies, we 
obtain constraints on parameters of the non-standard model. First, we study the model with positive 
tension of the brane without the 4D curvature term in the action. 

In the standard cosmology, in the radiation-dominated period, the radiation energy density ex- 
pressed as a function of t and of temperature of the Universe T is 

p r (t)=p r0 r(,), 4 = i ^- F , (32) 

Pr(T) = ^9*(T)T\ 
In the non-standard model, in the period of late cosmology we have 

I (t] „ Pror(z) p 2 

PrW ~ 2(H ty[(l + q )pi-3 + A(r(z)-l)n r /n m Y { ) 

~Pr(f) = ^*(T)T 4 

where by tilde we distinguish the non-standard case. Here z = z(t). 

The freezing temperature Tp of the reaction n <-> p is estimated as a temperature at which the 
Hubble parameter H is of order of the reaction rate G 2 F T F [17]. Time dependencies of the Hubble 



parameter in both the standard and non-standard cosmologies are the same H = l/2t. Using (32) 
and (33), we obtain the ratio of freezing temperatures in the standard and non-standard cosmologies 



fp = ( g.(f F ) [(1 + q )p 2 - 3 + Mrjz) - l)fl r /fl m ] * 
T F \g*(T F ) 2p 2 r(z)n r 

The mass fraction of A He produced in nucleosynthesis of the total baryon mass is 



(34) 



(n/p) f + 1 

where the subscript "f indicates that the ratio is taken at the end of primordial nucleosynthesis, 
(n/p)f ~ 1/7 [17]. The equilibrium value of the neutron-proton ratio (ti/p)t f = exp [—(m n — m p )/T F ] 
is very sensitive to the value of T F . Under variation of freezing temperature variation of X 4 is 

2 5T F 

6X4 " wm^ {n/v)Mv/n)F ^- (36) 

We constrain parameters of the model requiring that the difference between the calculated values of 
X 4 in the standard and non-standard models is within the experimental errors. The estimate of X 4 
obtained in the standard cosmology X 4 ~ 0.25 fits well the experimental value X 4 = 0.25 ±0.01 [18]. 
Using (36) we estimate the variation of the freezing temperature retaining X 4 within the experimental 
errors 

AT 1 

— - < 0.0255. 
T F 

Comparing T F and T F , we have 

1/6 



5T F _ ( g*(T F ) [(1 + q )p 2 - 3 + 4(r(z) - l)Q r /Q T 



T F \g*(T F ) 2p\{z)Q r 
From this relation follows an estimate 

g*(f F ) [(1 + q )p 2 - 3 + 4(r(z) - l)fi r /ft 



- 1. 



g*(T F ) 2p 2 r(z)n r 
where e = 0.164. It follows that 



< e (37) 



(l + q )p 2 -3 = O(n r ). (38) 
Introducing dimensionless ratio of scales fiM 2 jM 3 , we can write 

2 „ 2M3 

Substituting this relation in (37), we obtain an estimate 

A 2(l + g0 )(r-l) \ ml 2(1 + g ) 4rQ r / _ 2(1 + g )(r - 1) \ 

I £ m m r J M3 3fi m 3Q m I + 3n m r J' { ' 
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Now we can verify that in the non-standard model BBN takes place at the period of late cosmology. 
From condition of late cosmology (20), the transition point to late cosmology, can be bounded as 



z 4 < 



(40) 



In the most stringent case, for p ~ 10 _12 Gey, this yields z ~ 10 16 which is much larger than 
zbbn ~ 10 9H ~ 10 . In the radiation-dominated period, setting r(z) = 1, we obtain solution of the 
Friedmann equation as 



_ 4 ^2[(l + g y-3] W)2+ 4Q 



z ~ 



p 2 jiVL m Hot ' 



(41) 



From (41) and (40) it follows that the transition time to late cosmology is t ~ l/p ~ l0 12 GeV~ 1 . 
This is much smaller than characteristic time of nucleosynthesis 1 10 2 s, or lO^^GeV^ 1 . 

Let us turn to the model with 4D curvature term included in the action. Following the same steps 
as in the model with r c = 0, we find that production of Helium is within the experimental bounds if 



p\pr c + 1)(1 + g ) - 3 + 4(r - l)Q r /Q T 



2p 2 (fir c + l)rn r 
where e ~ 0.164. From (42) follows the constraint on p 2 



1 



< e, 



(42) 



3fi 



2(1 + go) 



<p 2 (l + /ir c )(l + g ) - 3 < 



ArQ r 



l + e) 



3fi 

r, 



r — 1 



2(1 + go) 



which can be written as 



p 2 (l + pr c ) = 



1 + go 



[i + o(n r )) 



(43) 
(44) 



Using the inequalities (43) we obtain the bounds on the (normalized) density of dark radiation p w0 



^ (1+ , rJ ( 1 _ E _?(I±<*> 

P y oi L m 



\ < P^ < ^^ 1+ )r f 1+e 

) fi p 2 v 



2(1 + g ; 

30 



(45) 



The combination 2(l + go)/3f2 m is of order unity, and within the existing uncertainties of cosmological 
parameters the sign of p w o is ambiguous. From (45) it follows that p r o/p > p w o/ '/-*■ At times at the 
beginning of nucleosythesis the effective number of degrees of freedom is q*(Tbbn) = 10.75, and we 
obtain r ~ 3.2. 

Solving the Friedmann equation (15) with r(z) = 1 in the radiation-dominated period of late 
cosmology, we have 



p 2 (pr c + 1) 

In (20) the period of late cosmology was defined as 



(46) 



z 4 < z 4 ~ 



p 



n r H 2 (l + pr c )' 



For the characteristic times of nucleosynthesis, tsBN ~ 10 24:26 GeV 1 , we have zbbn 
Requiring that zbbn *C z, we obtain 



10 



9-rlO 



Hi 
p 2 



Q z 4 

lL r 6 BBN 



(47) 
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In the most stringent case, for /i ~ 10 12 GeV, from (47) it follows that 

lir c < 10 25 , 

which is stronger than the bound (17). 

4.2 Model with negative tension of the brane 

First, we consider the model without the 4D curvature tern in the action. In the radiation-dominated 
period, solution of Friedmann equation is 

2 -^ 2(3 + (i +g0 y) 

P 2 P 2 M^m 

where we neglected the slow factor r(z). At the characteristic times of nucleosynthesis the first term 
in (48) is dominant, and we can neglect the second one. 

Following the same steps as in the model with positive tension of the brane, we obtain the ratio 
of the freezing temperatures in the standard and non-standard cosmologies 



T L= l g*{T F ) [3 + (l + q )p 2 ^ 

T F U.(7» 2 P 2 n r 



(49) 



Because (3 + (1 + q )p 2 )/2p 2 Q r is a large number, the ratio T F /T F is not close to unity, and the 
abundance of He calculated in this model is outside the experimental bounds. With the minimal 
value of the factor [3 + (1 + qo)p 2 ]/2p 2 Q r ~ 10 4 we obtain T F /T F = 4.64 and (n/p)f = 0.68. 

Using (48) and (32) we find that the times t and t the Universe has the same temperature T both 
in the standard and non-standard models are connected as 



t = t 



n rP 2 n 1/2 



p 2 (l + q ) + 3 



It follows that i < 10 2 t. In the non-standard model the interval of times at which the Universe is 
cooled from the temperature T ~ lMeV to O.lMeV between the beginning and the end of the BBN 
is ~ Is. The decay of neutrons during this time interval is negligible as compared with the standard 
model. The fraction X 4 can be calculated at the freezing point 

2(n/p) F 

x * = tt~\ — rr - °- 81 ' 50 

[n/p) F + 1 

and is considerably larger than the observational value. 

In the model with 4D curvature term included in the 5D action we considered two cases. If 
1 — \cr\r c > 0, the model is similar to the case without 4D term. As it was discussed above, this 
case the model yields the value of X 4 mach larger than the observational bound. If \cr\r c — 1 > 0, 
solutions of the Friedmann equations are similar to those in the model with positive tension of the 
brane and in this case the value of X 4 is within the experimental bounds if parameters of the model 
are restricted as 

p 2 (fir c - 1)(1 + g ) - 3 + 4(r - l)Q r /Q r 



2p 2 (fir c — l)r£l r 



- 1 



< e. 
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5 Conclusions 



In this note we discussed a class of warped cosmological solutions of the Friedmann equation in one- 
brane 5D models with infinite extra dimension. We considered models with positive and negative 
tensions of the brane. We studied separately the models with and without 4D curvature term on the 
brane included in the 5D action. We found constraints on parameters of the models which follow from 
consistency of predictions of the models on abundance of 4 i?e produced in primordial nucleosynthesis 
with the cosmological data. 

For numerical estimates we used the present-time values of the Hubble parameter Hq, fractions 
of cold matter and radiation Q m , Q r , deceleration parameter go, and the mass fraction X4 of A He of 
the total barion mass produced in primordial nucleosynthesis. 

It was found that in the models with positive tension of the brane, both with and without 4D 
curvature term, production of i He can be consistent with the cosmological data, if parameters of the 
models satisfy the constraint 



/^M 2 2(1 + g ) 

M*(nr c + 1) 3ft m 

Because 3f2 m /2(l + g ) ~ 1, it follows that 



< 0(fi r ). 



/iM 2 



fir c + 1, 



M 3 

where an upper bound on /ir c is 

Hr c < 10 25 . 

In the models with negative tension of the brane the results depend on whether 4D curvature term is 
included in the action of the 5D model or not. The model without 4D curvature term in the action 
yield abundance of A He considerably higher than the observational data. 

The model with 4D curvature term included in the 5D action can meet the observational data 
provided /ir c > 1. From the estimate of the abundance of 4 He follows a constraint on parameters 

/iM 2 2(1 + g ) 



< CQ r , (51) 



M 3 (/ir c - 1) 3fi m 
where C ~ 10. 

Further restrictions on parameters of the model with negative tension of the brane can be ob- 
tained, if Mp t is larger than the coupling at the 4D term in the action, i.e. if M 2 Z > r c M 3 , or 
fiMpjM 3 > fj,r c . In this case 



^ < ttM^/M 3 < [pr e - 1) ( 2(1 3 +J o) + CSKj 
2(1 + ^/2(1+^^ 



or 

2( 

tjii L m \ Oil L m 

Because 2(1 + q )/3Q m ~ 1, it follows that fir c ^> 1 and hence fiM^/M 3 ^> 1. 

Generally, numerical values of cosmological parameters depend on a model in which the experi- 
mental data were processed. Recently fits of parameters were made in RSII-like models with positive 
tension of the brane [21, 22, 23, 24]. The main cosmological parameters obtained in these fits do not 
considerably differ from those of the standard model used in the present note. With these parameters 
qualitative conclusions of the present paper remain valid. 

m 
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